Frequency and temporal effects in linear optical quantum computing 
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Typically linear optical quantum computing (LOQC) models assume that all input photons are 
completely indistinguishable. In practice there will inevitably be non-idealities associated with the 
photons and the experimental setup which will introduce a degree of distinguishability between 
photons. We consider a non-deterministic optical controUed-NOT gate, a fundamental LOQC gate, 
and examine the effect of temporal and spectral distinguishability on its operation. We also consider 
the effect of utilizing non-ideal photon counters, which have finite bandwidth and time response. 
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I. INTRODUCTION 

It has been shown that efficient quantum computa- 
tion can be implemented using only linear optical ele- 
ments, single photon sources and photon counters Q. 
This scheme is referred to as linear optical quantum com- 
puting (LOQC) and has recently been subject to several 
in-principle experimental demonstrations |^ y, |j, |^ - 

In the LOQC model qubits are encoded as the pres- 
ence in one of two modes of a single photon (dual-rail 
logic). Thus a logical zero can be written |0)l = |01) 
where \ij) is a two mode state in which i is the occu- 
pation number of the first mode and and j is the occu- 
pation number of the second mode. Similarly a logical 
one is written = |10). Often these modes are taken 
to be the vertical and horizontal polarization modes of a 
single spatial mode. The non-linearities which are neces- 
sary for quantum computation are introduced via a con- 
ditioning or post-selection process whereby the outputs 
from the system are only accepted if certain states are 
detected at extra ancillary outputs. This post-selection 
process results in non-deterministic quantum gate oper- 
ation (i.e. gates have a certain probability of failing), 
however, through the use of a teleportation protocol, it 
is possible in-principle to boost the probability of success 
arbitrarily close to unity . 

It is important to understand how critically non-ideal 
factors affect basic gate operation. Previous studies have 
looked at the effects of beamsplitter parameters ](^ , de- 
tection (7(1 and ancilla j^i] efHciencies on various gate de- 
signs. 

A key ingredient of the post-selection process is the 
presence of non-classical interference due to photon indis- 
tinguishability. Thus one would expect gate operation to 
be compromised if the input photons are distinguishable 
in some way. In this paper we consider such a problem by 
analyzing the effect of non-identical frequency distribu- 
tions and unsynchronized arrival times on the operation 
of a non-deterministic 2-qubit LOQC gate. Some previ- 



ous work has been done on the effect of unsynchronized 
arrival times in the context of quantum dot single pho- 
ton sources We begin by introducing a single pho- 
ton time/frequency representation. We then illustrate 
our techniques and some basic physical effects by the ex- 
ample of 2 photon interference at a 50/50 beamsplitter. 
In section IIVI we reapply these techniques to the more 
complicated case of the 2-qubit gate and also give con- 
sideration to frequency and temporal effects within the 
photon-counters. In section we conclude. 



II. REPRESENTATION OF SINGLE PHOTON 
STATES 

Typically when we employ single photon input states in 
the LOQC model we implicitly assume that all photons 
input into the system are completely indistinguishable. 
Hence, although all of the photons necessarily exhibit 
distributions in frequency space, it suffices to treat them 
all identically and express them in Dirac notation simply 
as |1), which ignores the exact form of their distribution. 
If, however, the input states are not completely indis- 
tinguishable then such a notation is insufficient and we 
must revert to a more general representation for single 
photon states which allows for their individual distribu- 
tions. This can be modeled by adopting the notation 



|1) 



\l)..dLU 



au^al duj 



(1) 



where a^j is a normalized function describing the dis- 
tribution of the photon in frequency space. is the 
single frequency, single-photon state. Similarly, aj^ is 
the single frequency creation operator. The integral is 
over all frequency space. In choosing this representa- 
tion of frequency distributed states, we have implicitly 
assumed that inputs are pure states. Consequently pulses 
are transform limited and the corresponding time-domain 
representation of any distribution can be determined by 
taking the inverse Fourier transform of the distribution 
function. Wc will discuss mixed states at the end of this 
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By choosing different forms for a^^, photons produced 
by various physical systems can be modeled. We now 
consider some examples. The output from an optical cav- 
ity is characterized by a Lorentzian distribution, which 
has the form 



1 



(2) 



where k is a parameter related to the bandwidth of the 
distribution. Thus we would expect that a single photon 
source, comprising a single photon emitter in an opti- 
cal cavity, would be approximately characterized by a 
Lorentzian distribution. We say approximately because 
we assume that there is no uncertainty in the time of 
emission of the photon by the photon emitter. In re- 
ality there would be some uncertainty associated with 
this and the distribution of the photon emitted from the 
cavity would be represented as a convolution of this dis- 



tribution and the cavity's response. However, assuming 
that the time-frame in which photons are placed into the 
cavity by the photon emitter is much smaller than the de- 
cay time of the cavity, a Lorentzian distribution is a valid 
approximation of the res pon se. Examples are quantum 
dot 01 and cavity QED based photon sources. The 
Lorentzian distribution has the property that, although 
two-sided in the frequency-domain, it is one-sided in the 
time-domain. 

Another possible single photon source is the non- 
degenerate parametric down converter, frequently used 
in LOQC experiments. This is a device which takes a sin- 
gle photon as an input and produces two photons, each 
with half the frequency of the original photon, but in 
orthogonal spatial or polarization modes, at its outputs. 
Considerable work has been done on the frequency dis- 
tributions of down-conversion sources [T^ Il3j | . The form 
of the output from a down-converter can be written 



|^o«t) - |0)a|0)f,- 



2x« 



■(|0^1 |UO_^>fc + |i^o 



(3) 



where x is the conversion efficiency of the down- 
converter, and K is again related to the bandwidth of the 
output photons. We assume that conversion is very weak 
such that X ^ K. This justifies dropping higher order 
terms involving higher photon numbers in Eq. The 
modes a and b correspond to the two spatial modes which 
form the outputs. It is evident from the form of Eq. 
that the down-conversion process does not always pro- 
duce photon pairs, since it is possible that both outputs 
will be in the vacuum state. In fact, with the assump- 
tion that X is small compared to k, the down-converter 
will produce vacuum outputs most of the time. However, 
because the down-converter produces photon pairs if it 
produces anything at all, a reliable photon source can be 
constructed by conditioning on the detection of a photon 
at one output and using the other output as the source. 
Having performed conditioning, the normalized form of 
the photon produced at the other down-converter output 
is 

- / ^ 2xK 

where it has been assumed that the intrinsic response 
time of the counter, r, obeys 1/t ^ k (see section llV|l . 
Unlike the Lorentzian distribution, this distribution is 
a symmetric, real function in both the time and fre- 



quency domains. Passing a photon of this form through 
an appropriate filter can produce a Gaussian distributed 
packet. The normalized form of such a distribution is 




where again k is related to the bandwidth of the distri- 
bution. 

In addition to modeling arbitrary frequency distribu- 
tions, we can also introduce an arbitrary time-delay onto 
a photon wave-packet. Using Fourier transform rules this 
can be modeled by introducing a complex exponential 
rotation factor into the frequency-domain distribution 
function. This corresponds to making the substitution 
Qfij e^^^a^;, where r is the time-shift parameter. Thus, 
as an example, a Gaussian distribution with a time-shift 
would be expressed as 

ac. = e^"^'H-e-^ (6) 

V K TT 

Finally, mixed states can be represented by introduc- 
ing a classical distribution of pure states. For example 
a single photon with inhomogeneous broadening of its 
frequency distribution can be represented by the density 
operator 



3 



dfle ~r / / duJduJ'a^+Qal,_^_^-l\l)^{l\^, (7) 



In the following analysis we will restrict ourselves to pure 
input states. 



III. THE BEAMSPLITTER 

To illustrate our method we begin by considering the 
operation of the beamsplitter, which is a fundamental 
building block in LOQC. We adopt the convention of the 
phase-asymmetric beamsplitter, shown in Fig. ^ which 
has the property that reflection off either surface or trans- 
mission from the 'black' surface result in no sign change, 
whereas transmission from the 'gray' surface results in a 
sign change. 




FIG. 1: The phase asymmetric beamsplitter. Reflection off 
either surface or transmission from the 'black' surface result in 
no sign change. Transmission from the 'gray' surface results 
in a sign change. 

The evolution of states through the beamsplitter is de- 
scribed in the Heisenberg picture by the equations of mo- 
tion 

Siout = V^a + ^1 -r/h 

bout = - Vl - (8) 

where a and b are the two spatial modes and r] is the 
beamsplitter's reflectivity. 



A. Calculating the beamsplitter output state 

We now describe how the Schrodinger picture evolu- 
tion of the state may be obtained from Eq. 0. First 
we express our input in terms of creation operators. For 
example, the input state with one photon at each beam- 
splitter input would be expressed as 

IV) = |l>a|l>b^atbV)a|0)b (9) 



To determine the corresponding output state for this in- 
put state we apply the time-evolution operator U, to the 
input state. We also insert the identity operator UW, 
after each creation operator. Finally we recognize that 
when both inputs are in the vacuum state, both outputs 
must also be in the vacuum state. Hence we can apply 
the simplification C/|0)a|0)b = |0)a|0)b. Thus the output 
state can be represented as 

\lp)out = UWin 

= ?7atbV)a|0)6 

= uk\u^u)h\u^u)\o)a\o)t 
= {uk^W){uh^u^)u\o)a\o)b 

= iUeL^U^)iUh^U^)\0)a\0)b (10) 

The problem is now reduced to one of finding the reverse 

time evolved creation operators U a' and Uh . We 
can determine these by first rewriting the equations of 
motion from Eq. (|SJ) as 

SLout = U''aU = + ^l-Tjh 

hout = U^hU = V^b - v/l-7?a (11) 

By applying U from the left and from the right of this 
expression we obtain 

a = ^UaU^ + v/1 - r]?7bC/t 

b = ^yC/bfjt - y/l - riUaU'' (12) 

This results in a linear system of equations in terms of the 
reverse time evolved annihilation operators, which can be 
inverted to determine tJkW and Uhtj\ and conjugated 

to determine J/a^C/^ and t/b^C/t. In this example the 
solution is 

J/a^C/t = V^a^ - ^1 - rih^ 

?7b^C/t = VT^at + (13) 

Substituting this result into the expression for the output 
state in Eq. (|10(l and setting rj = 0.5 results in 

\i^out) = (iat'-^b^')|0)a|0)b 

= -L|2)j0>b--^|0)j2)b (14) 

which displays quantum interference through complete 
suppression of the "|l)a|l)h" output term. This interfer- 
ence critically depends on the indistinguishability of the 
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photons. We now generaHze this calculation by allowing 
for differently distributed input states. We do this by 
making the substitution presented in Eq. Thus the 
input state is now expressed as 



l>a|l)6 



a,., a 



idc.)( / pSldLo)\0)a\0)b (15) 



It should be emphasized that this representation only 
allows for distinguishability in the time/frequency degree 
of freedom and does not allow for spectral entanglement. 
Upon applying the same procedure as before to this input 
state, we determine that the corresponding output state 
for the 50/50 beamsplitter is 



a^f3^,{\l)^\l)^,)a\0)bdujduj' 
|l)„'^f, dojduj' 

a^M^)a{\l)u\l)u')bduodu' (16) 



This result is completely general and allows us to substi- 
tute arbitrary frequency distribution functions. Without 
actually performing any substitutions however, we can 
gain much insight into the beamsplitter's behavior sim- 
ply by inspecting the form of this result. The expression 
contains four terms, the first and last of which correspond 
to the terms in the ideal situation (see Eq. ^5) where 
there are two photons present at one output and none 
at the other. The center two terms arise from photon 
distinguishability in which there is one photon present 
at each output. When and (3^ are equal the center 
two terms cancel leaving only the ideal terms. This is 
the expected result, since so long as the inputs are in- 
distinguishable ideal quantum behavior should ensue. In 
the opposing limit, when there is no overlap between a^^ 
and (3uj, no cancellation takes place and there will be a 
50% probability of measuring exactly one photon at each 
beamsplitter output. The photons then exhibit classical 
statistics. 

From this expression for the output state we can cal- 
culate the coincidence rate {NaNb), between the outputs 

{NaNb) = \ff |ac.n/3^fda;dc^' 



a^l3^*a^* 13^1 dujdu)' (17) 



For and (3^^ normalized we obtain 

{NaNb) = ^ ^ ^ // C(u,Puj*(^uj'*Pu,' dtuduj' 



Upon inspection we see that when and /3(^ are equal 
the two terms cancel, resulting in (NaNb) = as expected 
for indistinguishable photons. On the other hand, when 
there is no overlap between a^^ and the second term 
reduces to 0, since this term calculates the overlap be- 
tween the two distributions. Hence, assuming that the 
inputs are normalized, (NaNb) — 0.5 as expected for a 
classical distribution. 



B. Beamsplitter examples 

We now consider several examples of photons with dif- 
ferent frequency distributions and arrival times. 



1. Gaussian distributions 

Substituting Gaussian distributions (see Eq. jS]) into 
the general expression for (NaNb) obtained in Eq. H18(l 
results in 



1 



(NaNb) = - 



EK 



(19) 



where e and k are the bandwidth parameters associated 
with the two inputs. It is evident from this expression 
that when e — k, (NaNb) = 0. Alternately, if either e 
or K approach zero (but not both simultaneously) then 
{NaNb) = 0.5. This result is expected since a zero band- 
width distribution corresponds to a perfectly monochro- 
matic photon, which will be completely distinguishable 
from any non-monochromatic photon. Fig. [21 shows a 
plot of this behavior. We now consider the case where 




N.N:. 



(18) 



1.0 1.0 



FIG. 2: {NaNt) against the bandwidths of the two photons 
input to a beamsplitter, where both inputs are Gaussian dis- 
tributed and arrive simultaneously. The units for e and k are 
arbitrary, provided they are consistent. 

a time-shift r, is introduced into one of the beamsplitter 
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inputs. In this case we find that 



{Nam: 



1 



(20) 



where we have let £ = k for simphcity. Plotting this re- 
sult yields Fig. |21 which some readers will recognize as 
the famihar Hong-Ou-Mandel (HOM) ^ dip. This re- 
sult indicates that when no time-shift is present between 
identically distributed inputs, perfect quantum behavior 
arises. As the time-shift is increased and the two pho- 
tons become increasingly distinguishable, the expecta- 
tion value asymptotically approaches the classical result. 





-TIC 



FIG. 3: {NaNb) against the time-shift introduced into one of 
the inputs to a beamsplitter, where both inputs are Gaussian 
distributed and share the same bandwidth k. 



FIG. 4: (NaNt) against the time-shift introduced into one of 
the inputs to a beamsplitter, where both inputs are Lorentzian 
distributed and share the same bandwidth k. 



3. Down-conversion 

We now calculate (NaNi,) when using inputs derived 
from conditioned down-converters. The result is of the 
form 



(NaNb) 



1 



2eK 



(23) 



2. Lorentzian distributions 



where we have let x — 1 for simplicity. This result is 
identical to the Lorentzian case. Upon introducing a 
time-shift into one of the inputs we obtain 

(7V„iV,) = i-ie-2''l-l(l + «|r|)^ (24) 

where e = k. The behavior of this result is indistinguish- 
able from the Gaussian case, and, unlike the Lorentzian 
case, we again observe a smooth rather than a sharp dip. 



When using Lorentzian distributed inputs (see Eq. 
with different bandwidths, similar behavior arises to the 
Gaussian case. (NaNb) now has the form 



- 2 - (7W 



(21) 



This result behaves almost identically to the Gaussian 
case and has the same properties in the limits. 

Introducing a time-shift into one of the inputs results 



{NaNb] 



= -2k|t| 



(22) 



where we have again made the simplification e = k. Plot- 
ting this yields Fig. 0] which behaves similarly to the 
Gaussian case with the notable difference that the dip 
is no longer smooth, but rather sharp. This sharpness 
can be attributed to the discontinuity of the Lorentzian 
distribution in the time-domain. 



IV. THE CONTROLLED-NOT GATE 

A number of non-deterministic, heralded two qubit 
gates have now been proposed P, 0, 0|. We con- 
sider here the controUed-NOT (CNOT) gate proposed by 
Ralph et. al. This gate is simple in design and has 
previously been shown to exhibit good resilience to an- 
cilla inefficiencies when compared to several other LOQC 
CNOT gate implementations .Sj. 

The CNOT gate is the quantum equivalent of the clas- 
sical XOR-gate. The gate operates on two input qubits, 
the control (c) and target {t) , and implements the follow- 
ing transformations on the logical basis states \H) and 
\V) 



|^),„ = \H),\H)t 

W^n = \H)c\V)t 
W^n = \V)c\H)t 
W^n = \V)c\V)t 



mout = \H)c\H)t 
Want = \H)c\V)t 
Wont = |^^>c|^>t 
Wont = \V)c\H)t (25) 
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where \H) = \1)h\0)v and \V) = \0)h\1)v- Of course, 
the gate wih also operate on superposition states. 

The gate, shown in Fig. is implemented using eight 
beamsplitters, with reflectivities as indicated in the di- 
agram {r]i — yys). In addition to the control and target 
inputs the gate employs two vacuum [vi, V2) and two 
ancilla (oi, 02) inputs. The gate is non-deterministic, 
in that the success of the gate is conditional upon de- 
tecting no photon or exactly one photon at some of the 
outputs, labeled "0" and "1" respectively. Ideally the 
success probability is approximately 5%. 



A. Modeling the CNOT gate 

In our study of the CNOT gate we will only consider 
photon distinguishability effects and ignore other possible 
experimental errors. Using the same approach as for the 
beamsplitter we derive a general expression for the input 
into the CNOT gate, where the inputs are subject to 
frequency distributions. The form of the input state is 



rT 



7( / «^c|,_c^ dw^ / (T^t^^^dtj^ + 



(T^iy ,^ duj 



(26) 



duj ! 



where a, f3, 7 and S determine the input superposi- 
tion, and CTtj and e^: are the frequency distribution 
functions of the control, target and ancilla inputs respec- 
tively. 

Using this expression for a completely general input 
state, we wish to determine the corresponding output 
state. We use a similar approach as previously for the 
beamsplitter. Firstly, we derive the equations of motion 
for the system by recursively using the equations of mo- 
tion for a single beamsplitter. This system of equations 
is then inverted to determine the reverse time evolved 
creation operators for each of the input modes. These 
operators are used to evolve the input state given in Eq. 
(I26|l . We omit the resulting expression due to its com- 
plexity. 



B. Modeling photon counters 

Unlike the beamsplitter we considered previously, the 
CNOT gate involves a conditioning process which must 
also be modeled. This is performed using photon coun- 
ters. In a standard model, in which frequency distribu- 
tions of inputs are ignored, the conditioning process could 



be modeled by applying a series of projectors represent- 
ing the desired measurements. In the case of the CNOT 
gate this would be expressed as 



i\0U{0Um\0)vA0Um) 



(11 



II) 



(11 



(27) 



This sequence of projectors effectively discards terms 
which contain photons at any of the vacuum outputs or 
not exactly one photon at either of the ancilla outputs, 
the desired post-selection procedure. Clearly this se- 
quence of projectors is incompatible with our frequency- 
distributed model and must be modified. 

In the frequency distribution picture we modify the 
projector from Eq. (|27|) to 



(lO).i(OkJ® (|0).,(0|.J® ( / \l)a^AMa^,udiU 



(28) 



This projector describes post-selection using ideal photon 
counters. We now generalize the model further to accom- 
modate for two important real-life effects, the counters' 
frequency and time resolution. 

Any photon counter will invariably be sensitive to only 
a finite range of frequencies. More generally we can intro- 
duce a function, say ^uj, which describes this frequency 
response. We model the counter's frequency selectivity 
with the introduction of a filter prior to the projector. 
The filter is described by the equations of motion 



(29) 



In this model the counter's input enters through a while 
b is kept in the vacuum state. Uout then contains the 
transmitted (or accepted), and bout the reflected com- 
ponent of the input state. The function varies be- 
tween and 1, where indicates no sensitivity, and 1 
complete sensitivity, to the respective frequency compo- 
nent. We would expect that this function would gener- 
ally have a Gaussian-like distribution. The components 
which are reflected out of the 6-mode are 'unseen' by the 
observer. Hence this filtering process results in mixing ef- 
fects. Therefore, following this filtering stage the state of 
the system must be expressed in density operator form, 
where the reflected modes are traced out to produce a re- 
duced density operator describing the observed system. 
Thus 



PfiUered — tr^ ( IV'/iite 



ed){tpfiltered\) (30) 

]) is then applied to 



The integral form projector (Eq. 
this reduced density operator. 

In addition to flnite frequency selectivity, in any ex- 
periment the photon counters will only be 'open' for a 
finite time period. We model this time selectivity in an 
analogous manner to our modeling of frequency selectiv- 
ity, except that the equations of motion describing the 



FIG. 5: The layout of the simphfied CNOT gate in the LOQC model 



filtering process are now in the time-domain rather than 
the frequency-domain 

aout,T = y/rh^T + - rirhr 

hout,T = y/rp^T - - Vrar (31) 

where rjr is a time dependent function, varying between 
and 1, which describes the time-response of the counter. 
In general we would expect this function to be approx- 
imately rectangular, mimicking the opening and closing 
of the shutter. 



C. Characterizing CNOT gate behavior 

Having developed an expression for the output state 
of the CNOT gate for a general input state and photon 
counter response we wish to characterize the behavior of 
the gate as various parameters are changed. We use two 
quantifiers to do this, the fidelity and success probability 
of the gate. 

The fidelity is a measure of how close the actual output 
state and expected output state of the gate are, where the 
expected output state is defined according to the logical 
transformations given in Eq. H25|) . and we assume that 
all inputs share identical distributions and arrive simul- 
taneously. Hence, the fidelity quantifies how accurately 
the gate is performing the desired logical transformation. 
Mathematically this is defined as 

r^ (fpexplPoutl^exp) /'oo^ 

tr[pout) 

where \4>)exp is the normalized expected output state and 
Pout the density operator of the actual output state. The 
trace is introduced into the denominator to normalize 
Pout, thereby decoupling F from the success probability. 
Due to normalization we expect F to vary between and 



1, where 1 corresponds to ideal gate behavior, and to 
completely non-ideal behavior. We expect to see F = 1 
when all of the input photons are indistinguishable and 
the counters are behaving ideally (i.e. they have infinite 
bandwidths and time windows). 

We should caution that the use of a 1-dimensional mea- 
sure such as fidelity to characterize our gate rolls a num- 
ber of effects into a single number. In particular Itpexp) 
is defined as the expected output state when the input 
states have identical distributions and arrive simultane- 
ously. This has the effect that if, for example, a time-shift 
were introduced into one of the inputs then F would suf- 
fer as a result of two separate effects. The obvious effect 
is that the gate would no longer be operating ideally and 
thus the transformation implemented by the gate would 
be different to what is expected. The less obvious effect 
is that even if the logical transformation implemented 
by the gate were not changed, F would still suffer since 
\''Pexp) is defined as the output state when no time-shifts 
are present. Hence we are calculating the overlap of a 
state which is displaced temporally with one which is not. 
There is no obvious way to decouple the measurement of 
these two effects. On the other hand we should remem- 
ber that ultimately multiple quantum gates are intended 
to be integrated into quantum circuits. In the context of 
a circuit whereby the output from our gate is feeding into 
the input of another, both of the aforementioned effects 
are degrading and it is therefore valid that F takes both 
into account. 

The success probability measures how often the condi- 
tioning process, and hence gate operation, succeeds. This 
is defined as 

P = tiipout) (33) 

Thus the success probability is simply a reflection of the 
normalization factor of the output state. Recall that in 
the ideal case the CNOT gate succeeds approximately 
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5% of the time. Hence in the ideal hmit we expect P to 
approach 0.05. 

Using these two measures we can effectively character- 
ize how the behavior of the CNOT gate varies for any 
given input superposition as the input distribution func- 
tions and counter parameters are varied. Gate perfor- 
mance is input state dependent, so instead of considering 
a particular superposition we consider the worst-case sce- 
nario whereby we search across all possible input super- 
positions and determine the minimum values for both the 
fidelity and success probability (Fmin and Pmin)- Hence 

i^mm = niin(F) V a,/3,7,(5 
P™„=min(P) V a, 13,-1,5 (34) 

In our studies we perform the searches using a Monte 
Carlo approach. 



D. CNOT gate with time-shifted Gaussian inputs 

From the beamsplitter examples in section IHII we ob- 
serve that qualitatively the results are similar across the 
different input distributions considered. For simplicity 
we now restrict ourselves to considering Gaussian dis- 
tributed inputs and postulate that different distributions 
will yield similar results. 

Time-shifted Gaussian distributions (see Eq. [S]) were 
substituted for all of the frequency distribution functions 
in the expression for the CNOT gate's input state, with 
all bandwidth parameters equal. Hence, the inputs into 
the gate differed only in their arrival time. We sub- 
stituted flat distributions for the photon counters' fre- 
quency and time response functions. Thus the counters 
were ideal. Using these parameters, Fmin and Pmin were 
calculated as each of the time-shift parameters were var- 
ied. Fig. El illustrates the results when the Tcontroi (top) 
and Tanciiia (bottom) parameters are varied. It should 
be noted that, in the context of a search across all possi- 
ble input states, varying Ttarget is equivalent to varying 
^control, since beamsplitters Bl and B4 are invisible to 
such a search and thus the gate becomes completely sym- 
metric. Similarly, varying Tcontroi and Ttarget in unison is 
equivalent to varying Tanaiia- 

It is immediately obvious that the fidelity plots behave 
as expected in the limits. Specifically, when no time-shift 
is present in the input in question, perfect fidelity is ob- 
served. As the time-shift in the input approaches infinity, 
the fidelity asymptotically approaches zero {i.e. the out- 
put of the gate no longer bears any resemblance to what is 
expected). Similarly, when no time-shift is present we ob- 
serve the ideal-case success probability of approximately 
5%, which decreases as the time-shift is increased. Note 
however, that in general the success probability does not 
asymptotically approach zero. In fact, when either the 
control or target inputs are displaced (but not both si- 
multaneously) Pmin is completely unaffected. If P is 
plotted for specific superpositions we often see that P 



actually increases as time-shift parameters are increased. 
Current estimates place the fault-tolerant threshold for 
LOQC at the 1% level Our results imply that time 
synchronization needs to be to better than 1/10 of the 
inverse frequency bandwidth (1/k) to satisfy this thresh- 
old. Equivalently we can estimate that inhomogeneous 
temporal broadening (time jitter) must have a standard 
deviation less than this figure to reach fault tolerance. 
This result is somewhat better in terms of tolerable iit- 
ter than that obtained b y K iraz et. al. 9] for Lorentzian 
inputs to the Knill gate |15j . 

Although we have performed our calculations explic- 
itly as time-displaced inputs, transform symmetry means 
completely equivalent graphs would be generated as a 
function of frequency displacements and so we can simi- 
larly estimate that inhomogeneous frequency broadening 
should be less than 1/10 of the inverse temporal band- 
width. More generally any variation in the distinguisha- 
bility of the input qubits can be mapped to the results 
presented here. This is because it is inherently the degree 
of distinguishability between input states which deter- 
mines the effectiveness of the gate's operation, not the 
specific nature or source of the distinguishability. The 
conclusion is that although the results presented are very 
specific, their interpretation is very general and the pa- 
rameters Tcontroi, Ttarget and Tanciiia Can bc morc broadly 
interpreted as 'distinguishability factors' than specifically 
as time-shifts. We note that although the results pre- 
sented are for pure inputs states, results for mixed input 
states (see Eq. |7]) could be obtained by performing a 
weighted average over the results shown. 



E. CNOT gate with non-ideal photon counters 

We finally consider the effect upon the operation of 
the CNOT gate of utilizing photon counters which have 
finite frequency or time resolution. Initially we examine 
frequency resolution effects by substituting a Gaussian 
distribution for the frequency response function /i^, and 
leaving the time-range over which the counters operate 
as infinite. We assume that all of the input states are 

indistinguishable [i.e. Tcontroi = Ttarget = Tanciiia = 0). 

Fig. [7| shows a plot of the worst-case fidelity (a) and suc- 
cess probability (b) against photon counter bandwidth. 
As the counter bandwidth approaches infinity, both the 
fidelity and success probability approach the ideal limit. 
In the opposing limit, as the counters becomes completely 
selective, the success probability approaches zero and the 
fidelity approximately 0.2. In other words, the gate never 
succeeds and the gate's output is essentially random. 

Next we examine time resolution effects by assum- 
ing that the photon counters have infinite bandwidth, 
but operate only over a finite time period. We model 
this by substituting a rectangular function, ranging from 
~Ty,indow to Twindow, for the time response function rir. 
This is intended to model the behavior of a counter with 
a shutter with very fast reaction time. Hence the tran- 
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FIG. 6: Plots of Fmin (a,b,d,e) and Pmin (c,f) for the CNOT gate as the time-shift parameters of the control (top) and ancilla 
(bottom) inputs are varied. Plots (b) and (e) are fidelity plots, zoomed in on the region of high fidelity {F > 0.85). All inputs 
are Gaussian distributed with bandwidth k. 




FIG. 7: Fmin (a) and Pmin (b) against the bandwidth of the photon counters. The counters operate over an infinite time-range. 
All inputs are Gaussian distributed with bandwidth k and simultaneous arrival time. 



sient behavior of the shutter is ignored. The resuhs are 
shown in Fig. |S1 which exhibits similar behavior to the 
plots against the counters' bandwidths. 

Legero et. al. recently performed a HOM-like ex- 
periment in which photons were intentionally introduced 
with a relative time-shift. It was found that if the time 
resolution of the photon counters was sufficiently narrow 
compared to the photons' length, photon distinguishabil- 
ity could be compensated for by 'zooming in' on the time 
period where the overlap between the photons' wave- 
functions was maximal and the photons were effectively 
indistinguishable. One might question whether a similar 
approach is possible in the context of the CNOT gate 
to improve fidelity at the expense of success probability. 
Fig. 1^1 illustrates this concept. Upon investigation it is 
evident that in general this approach does not work for 
the implementation of the CNOT gate described. This 
is due to photon- number ambiguity at the counters. In 
the HOM experiment it is known, due to conservation of 
the number of photons in the system, that if both coun- 
ters 'click' exactly once then exactly one photon must 
have been incident on each counter. In the context of 



the CNOT gate we cannot make such assumptions due 
to the extra ancillary photons which make it possible 
that anywhere between zero and three photons will be 
incident upon the counters. Hence, if a counter with a 
finite time window clicks once then this does not neces- 
sarily imply that exactly one photon was incident upon 
the counter. It is possible that two or even three pho- 
tons were incident, but were simply not seen due to the 
finite time window. For this reason we postulate that in 
general it is not possible to improve the operation [i.e. 
fidelity) of the CNOT gate by manipulating the time 
window of the conditioning photon counters. We note 
that this argument does not apply to current demonstra- 
tions of the CNOT gate which operate in coincidence, as 
photon-number ambiguity is no longer an issue. Calcula- 
tions on the O'Brien et. al. gate !^ 1^ show that fidelity 
can be improved arbitrarily close to unity, at the expense 
of success probability, for arbitrary temporal mismatch 
between the inputs. The reasoning applied here to tem- 
poral filtering, when applied in the frequency domain, 
leads to the well known advantages of frequency filtering 
in experiments which operate in coincidence. 
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FIG. 8: Fmin (a) and Pmin (b) against the size of the time- window of the photon counters. The counters have infinite bandwidth. 
All inputs are Gaussian distributed with bandwidth k and simultaneous arrival time. 



At 




FIG. 9: A plot of the arriving photon wave-packets, where 
there is a time displacement of the ancilla photons relative to 
the other two photons. The region At is where wave- function 
overlap is maximal and the photons are almost indistinguish- 
able. 



CONCLUSION 



looked at how non-classical interference at a beamsplitter 
is affected by differing frequency distributions and time 
shifts. As an example of an LOQC circuit we examined 
the effect of time-shifts on a non-deterministic, heralded 
CNOT gate. We concluded that fidelities higher than 
99% required photon synchronization to better than 10% 
of the inverse bandwidth. Similar results will apply to 
frequency shifts. Finally we examined the effects of non- 
ideal photon counters upon the operation of the CNOT 
gate, specifically considering the time and frequency re- 
sponse of the counters. We found that in general the 
fidelity of the gate diminishes as either the frequency or 
time bandwidth of the counters decreases. Furthermore, 
the gate's fidelity cannot be improved through manipula- 
tion of the photon counters' frequency or time-responses. 
This is in contrast to the situation with gates which oper- 
ate in coincidence where such manipulations can be quite 
effective. 
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